144216XZ 



U.S. Patent Application For 



NON-UNIFORM VIEW WEIGHTING 
TOMOSYNTHESIS METHOD AND 
APPARATUS 



By: 

Baojun Li 
Bernhard Erich Hermann Claus 
Gopal B. Avinash 
Stephen W. Metz 
Jiang Hsieh 



EXPRESS MAIL MAILING LABEL 

NUMBER: EL 982 235 490 US 

DA TE OF DEPOSIT: November 26, 2003 



Pursuant to 37 C.F.R. §1.10, I hereby certify that I am personally depositing this paper 
or fee with the U.S. Postal Service, "Express Mail Post Office to Addressee " service on 
the date indicated above in a sealed envelope (a) having the above-numbered Express 
Mail label and sufficient postage affixed, and (b) addressed to the Commissioner for 
Patents, Mail Stop Patent Application, P.O. Box 1450, Alexandria, VA 22313-1450. 



November 26, 2003 
Date 




144216XZ 



NON-UNIFORM VIEW WEIGHTING TOMOSYNTHESIS 
METHOD AND APPARATUS 

5 BACKGROUND OF THE INVENTION 

The present invention relates generally to tomosynthesis imaging systems. More 
particularly, the invention relates to a technique for correcting or avoiding certain 
artifacts and generally improving tomosynthesis images by use of non-uniform view 
weighting. 

10 

Tomosynthesis systems have developed over the past decades and now offer 
significant advantages for many types of imaging, such as in the medical diagnostics 
field. In general, X-ray tomosynthesis involves the use of low dose X-ray radiation to 
produce a series of images acquired over a range of X-ray beam orientations relative to 

15 an image object. In currently available systems the object is positioned in front of a 

digital detector and the X-ray source is moved to various positions to produce the series 
of images. A number of such images may be produced, typically in excess of 50 or 60. 
The detector collects electrical data representative of the depletion of a charge at 
individual pixel locations resulting from continuation of the X-ray radiation at those 

20 locations by intervening objects, such as the features of the images subject, a patient in 

the medical diagnostics context, for example. Acquiring images of the subject from a 
different orientations of the X-ray beam allows depth information to be incorporated 
into the final 3D image. The depth information is unavailable in conventional projection 
X-ray imaging, making tomosynthesis attractive for identifying specific features of 

25 interest and their general location within the subject. 

Amorphous silicon flat panel digital X-ray detectors are currently available for 
tomosynthesis imaging. In general, however, any X-ray detector that provides a digital 
projection image may be used. These may include, for example, charge coupled device 
30 (CCD) arrays, digitized film screens, or other digital detectors, such as direct conversion 

detectors. The low electronic noise and fast read-out times of such detectors enable 
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acquisitions with many projections at low overall patient dose as compared with 
competing detector technologies. 

Following acquisition of the image data, tomosynthesis techniques include 
5 reconstruction of images at various "slices" through the subject. Reconstruction 

algorithms permit reconstruction of many such slices at different spatial planes, typically 
parallel to the imaging plane of the detector. Such slices contain different anatomies 
located at various heights above the detector with underlying and overlying structures 
being generally suppressed. The generation of slices from projection images typically 
10 contains a chain of processing and operations. It should be noted that the term 

"projection" generally refers to a specific geometry or positioning of the X-ray source 
with respect to the subject and detector, many such projections being used in generating 
the slices as noted above. 

15 Processing and operations used to generate slices from projection images in 

tomosynthesis include, but are not limited to, standard pre-processing steps, special pre- 
processing steps, reconstruction steps, and post-processing operations. Standard pre- 
processing operations include detector corrections, such as for gain, offset, bad pixels in 
the detector, and so forth. These may also include correction for geometry distortions, 

20 log transformation to store a "film-like look" and so forth. Special pre-processing steps 

may include bad detector edge correction, padding, beam hardening correction, off-focal 
radiation correction, reference normalization, and so forth. Reconstruction steps and 
algorithms may also vary. Many such algorithms are based on filtered back-projection 
principles, such as shift and add techniques, generalized filtered back-projection 

25 techniques, order statistics back-projections, and so forth. Algorithms based on back- 

projection filtering principles in which back-projection is performed first followed by 
2D/3D filtering, and algorithms based on minimum-norm solutions are also available. 
These include algorithms known in the field as ART, DART, MITS, TACT, Fourier- 
based reconstruction, objective function-base reconstruction, ML, MAP, and so forth, 

30 and combinations of these. Post-processing may include various types of image and 
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contrast enhancement, such as tissue equalization, thickness compensation, brightness 
and white balancing, and additional artifact management routines. 

In general, back-projection is the favored process to generate tomosynthesis 
5 slices, and is typically used in many of the reconstruction algorithms listed above. In 

back-projection, pixel intensity values are assigned to computed pixels of slices at the 
various levels above the imaging plane to form the slice images from the projection data. 
This is generally done by assigning values of pixels at each slice by dividing the 
intensity of the pixels in the projection data by the number of slices to be reconstructed. 

10 This uniform view weighting, however, causes many problems and inconsistencies in 

the reconstructed images. These issues may, then, result in artifacts such as wavy 
patterns in the reconstructed images, general fall off of intensity values near edges of the 
images, and ghost-like images at elevations above, image plane where the subject cannot 
have been located (i.e. beyond the physical limits of the subject being imaged). There is 

15 a need, therefore, for improved techniques for tomosynthesis image data processing. 

There is a particular need for a technique which accounts for a wide range of factors 
which can render reconstructed images inconsistent or otherwise degrade the images. 
Needs exist, for example, for techniques that reduce the fall-off of data near the edges of 
slice images, that reduce the assignment of values beyond the limits of the subject, and 

20 that reduce the deformity or inconsistency of the resulting data that can result from many 

factors affecting the underlying projection data. 

BRIEF DESCRIPTION OF THE INVENTION 

25 The present invention provides novel techniques for processing tomosynthesis 

data designed to respond to these needs. These techniques may be used in a wide range 
of processing steps, but are particularly well-suited to back-projection processes of the 
type used in many tomosynthesis image reconstruction approaches. The techniques 
provide for weighting of views in tomosynthesis which is non- linear for one or several 

30 reasons. The technique may be based, for example, upon creation of weighting maps 

which may include weights that are specific for each pixel in the various slices. A count 
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map may be included in a processing, along with a slice weighting function, and a 
projection weighting function. These maps may be computed in advance of imaging 
and used in reconstruction, or may be adapted for specific imaging implementations and 
exams. 

5 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a diagrammatical representation of an exemplary tomosynthesis 
system wherein non-uniform view weighting is implemented in accordance with the 
10 present techniques; 

Fig. 2 is a diagrammatical representation of a physical implementation of the 
system of Fig. 1; 

15 Fig. 3 is an elevational diagram of radiation paths during a tomosythesis image 

acquisition series illustrating issues involved in weighting of slice data accounted for 
through the present technique; 

Fig. 4 is a plan view of an exemplary count map used to non-uniformly weight 
20 image data in accordance with the present technique; 

Fig. 5 is a graphical illustration of a slice weighting technique used to account 
for physical limits of a subject being imaged; and 

25 Fig. 6 is an illustration of projection weighting masks used in the present 

technique for accounting for variations in response from a tomosynthesis system in 
various projections. 
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DETAILED DESCRIPTION OF SPECIFIC EMBODIMENTS 

Fig. 1 is a diagrammatical representation of an exemplary tomosynthesis 
system, designated generally by the reference numeral 10, for acquiring, processing 
and displaying tomosynthesis images, including images of various slices through a 
5 subject of interest in accordance with the present techniques. In the embodiment 

illustrated in Fig. 1, tomosynthesis system 10 includes a source 12 of X-ray radiation 
which is movable generally in a plane, or in three dimensions. In the exemplary 
embodiment, the X-ray source 12 typically includes an X-ray tube and associated 
support and filtering components. 

10 

A stream of radiation 14 is emitted by source 12 and passes into a region of a 
subject, such as a human patient 18. A collimator 16 serves to define the size and 
shape of the X-ray beam 14 that emerges from the X-ray source toward the subject. A 
portion of the radiation 20 passes through and around the subject, and impacts a 
15 detector array, represented generally by reference numeral 22. Detector elements of 

the array produce electrical signals that represent the intensity of the incident X-ray 
beam. These signals are acquired and processed to reconstruct an image of the 
features within the subject. 

20 Source 12 is controlled by a system controller 24 which furnishes both power 

and control signals for tomosynthesis examination sequences, including position of 
the source 12 relative to the subject 18 and detector 22. Moreover, detector 22 is 
coupled to the system controller 24 which commands acquisition of the signals 
generated by the detector 22. The system controller 22 may also execute various 

25 signal processing and filtration functions, such as for initial adjustment of dynamic 

ranges, interleaving of digital image data, and so forth. In general, the system 
controller 24 commands operation of the imaging system to execute examination 
protocols and to process acquired data. In the present context, the system controller 
24 also includes signal processing circuitry, typically based upon a general purpose or 

30 application-specific digital computer, associated memory circuitry for storing 
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programs and routines executed by the computer, as well as configuration parameters 
and image data, interface circuits, and so forth. 

In the embodiment illustrated in Fig.l, the system controller 24 includes an X- 
5 ray controller 26 which regulates generation of X-rays by the source 12. In particular, 

the X-ray controller 26 is configured to provide power and timing signals to the X-ray 
source. A motor controller 28 serves to control movement of a positional subsystem 
32 that regulates the position and orientation of the source with respect to the subject 
and detector. The positional subsystem may also cause movement of the detector, or 

10 even the patient, rather than or in addition to the source. It should be noted that in 

certain configurations, the positional subsystem 32 may be eliminated, particularly 
where multiple addressable sources 12 are provided. In such configurations, 
projections may be attained through the triggering of different sources of X-ray 
radiation positioned accordingly. Finally, in the illustration of Fig. 1, detector 22 is 

15 coupled to a data acquisition system 30 that receives data collected by read-out 

electronics of the detector 22. The data acquisition system 30 typically receives 
sampled analog signals from the detector and converts the signals to digital signals for 
subsequent processing by a computer 34. Such conversion, and indeed any 
preprocessing, may actually be performed to some degree within the detector 

20 assembly itself. 



Computer 34 is typically coupled to the system controller 24. Data collected 
by the data acquisition system 30 is transmitted to the computer 34 and, moreover, to 
a memory device 36. Any suitable type of memory device, and indeed of a computer, 

25 may be adapted to the present technique, particularly processors and memory devices 

adapted to process and store large amounts of data produced by the system. 
Moreover, computer 34 is configured to receive commands and scanning parameters 
from an operator via an operator workstation 38, typically equipped with a keyboard, 
mouse, or other input devices. An operator may control the system via these devices, 

30 and launch examinations for acquiring image data. Moreover, computer 34 is adapted 

to perform reconstruction of the image data as discussed in greater detail below. 
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Where desired, other computers or workstations may perform some or all of the 
functions of the present technique, including post-processing of image data simply 
accessed from memory device 36 or another memory device at the imaging system 
location or remote from that location. 

5 

In the diagrammatical illustration of Fig. 1, a display 40 is coupled to the 
operator workstation 38 for viewing reconstructed images and for controlling imaging. 
Additionally, the image may also be printed or otherwise output in a hardcopy form 
via a printer 42. The operator workstation, and indeed the overall system may be 

10 coupled to large image data storage devices, such as a picture archiving and 

communication system (PACS) 44. The PACS 44 may be coupled to a remote client, 
as illustrated at reference numeral 46, such as for requesting and transmitting images 
and image data for remote viewing and processing as described herein. It should be 
further noted that the computer 34 and operator workstation 38 may be coupled to 

15 other output devices which may include standard or special-purpose computer 

monitors, computers and associated processing circuitry. One or more operator 
workstations 38 may be further linked in the system for outputting system parameters, 
requesting examinations, viewing images, and so forth. In general, displays, printers, 
workstations and similar devices supplied within the system may be local to the data 

20 acquisition components or, as described above, remote from these components, such 

as elsewhere within an institution or in an entirely different location, being linked to 
the imaging system by any suitable network, such as the Internet, virtual private 
networks, Ethernets, and so forth. 

25 Referring generally to Fig. 2, an exemplary implementation of a tomosynthesis 

imaging system of the type discussed with respective of Fig. 1 is illustrated. As 
shown in Fig. 2, an imaging scanner 50 generally permits interposition of a subject 18 
between the source 12 and detector 22. Although a space is shown between the 
subject and detector 22 in Fig. 2, in practice, the subject may be positioned directly 

30 before the imaging plane and detector. The detector may, moreover, vary in size and 

configuration. The X-ray source 12 is illustrated as being positioned at a source 
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location or position 52 for generating one or a series of projections. In general, the 
source is movable to permit multiple such projections to be attained in an imaging 
sequence. In the illustration of Fig. 2, a source plane 52 is defined by the array of 
positions available for source 12. The source plane 54 may, of course, be replaced by 
5 three-dimensional trajectories for a movable source. Alternatively, two-dimensional 

or three-dimensional layouts and configurations may be defined for multiple sources 
which may or may not be independently movable. 

In typical operation, X-ray source 12 projects an X-ray beam from its focal 
point toward detector 22. A portion of the beam 14 that traverses the subject 18, 
results in attenuated X-rays 20 which impact detector 22. This radiation is thus 
attenuated or absorbed by the internal features of the subject, such as internal 
anatomies in the case of medical imaging. The detector is formed by a plurality of 
detector elements generally corresponding to discrete picture elements or pixels in the 
resulting image data. The individual pixel electronics detect the intensity of the 
radiation impacting each pixel location and produce output signals representative of 
the radiation. In an exemplary embodiment, the detector consists of an array of 2048 
x 2048. Other detector configurations and resolutions are, of course, possible. Each 
detector element at each pixel location produces an analog signal representative of the 
impending radiation, that is converted to a digital value for processing. 

Source 12 is moved and triggered, or distributed sources are similarly 
triggered, to produce a plurality of projections or images from different source 
locations. These projections are produced at different view angles and the resulting 
25 data is collected by the imaging system. In an exemplary embodiment, the source 12 

is positioned approximately 180cm from the detector, in a total range of motion of the 
source between 31cm and 131cm, resulting in a 5° to 20° movement of the source 
from a center position. In a typical examination, many such projections may be 
acquired, such as in excess of 50 or 60, although this number may vary. 

30 
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Either directly at the imaging system, or in a post-processing system, data 
collected by the system is manipulated by performing mathematical operations 
designed to compute the locations at which objects are located in the subject that 
produced the attenuation of the X-ray beam sensed by the detector. This overall 
5 process is generally referred to as back-projection. It has been found that a number of 

problems exist in conventional back-projection as applied to tomosynthesis imaging. 
For example, such techniques do not typically account for the fact that different 
regions of different slices are traversed by radiation for certain projections and not 
others. Because the back-projection techniques typically assume that all reconstructed 

10 slices contribute equally to X-ray attenuation, however, degradation of image quality 

for certain slices occurs because this assumption is incorrect. Similarly, conventional 
techniques tend to result in an assumption that locations physically outside of a 
subject may have contributed to the X-ray attenuation. Similarly, the techniques 
assume that all projections, corresponding to various source and/or detector locations, 

15 are equally valid for all pixels in the projections. As discussed in greater detail below, 

these assumptions are corrected by the present technique by providing non-uniform 
weighting of views in the image reconstruction process. Details regarding the specific 
steps in image reconstruction are not provided herein separately, as these steps may be 
generally conventional in nature. 

20 

Fig. 3 illustrates an elevational view of a tomosynthesis system in which two 
projection configurations are illustrated. The projection configuration, generally 
represented by reference numeral 56, is illustrated as including two source positions 
58 and 60, with source position 60 generally corresponding to a centralized source 

25 position with respect to detector 22, and source position 58 being offset from the 

center position. Each source position produces a beam 62 generally in a fan shape, 
which in actual practice may be a cone or controlled cone. The beam fans 62 impact 
detector 22 generally at an imaging plane 64. Above the imaging plane 64, various 
slices 66 may be constructed through the subject 18 to provide an indication of the 

30 internal features of the subject. 
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As can be seen in Fig. 3, the radiation emanating from source positions 58 and 
60 intersect one another differently in different regions of the imaging volume. 
Because slices 66 divide the imaging volume, then, the radiation intersects different 
slices differently and in spatially different locations. As shown in Fig. 3, for example, 
5 in an intersecting region 68, radiation from source position 58 fully intersects 

radiation from source position 60. In an other region, generally a non-intersecting 
region 70, however, an upper slice 66 is not traversed by radiation originating at 
source position 58, but is traversed by radiation originating from source position 60. 

Conventional back-projection reconstruction techniques would assume that 
attenuation of radiation from the various slices designated for reconstruction receive a 
contribution from objects at each slice position. However, as can be seen in Fig. 3, 
this is not the case. That is, all slices illustrated in the figure do contribute to 
attenuation of radiation from source positions 58 and 60 in the intersecting region 68, 
although within region 70, only radiation from source position 60 is attenuated by the 
features visible in the slices. It may also be noted from Fig. 3 that the surface area or 
portion of each projection in which the radiation intersects and does not intersect is 
not equal. That is, each slice has different regions of intersecting and non-intersecting 
radiation. It should also be noted that, in practice, many more projections are 
generally generated, resulting in many such regions of intersecting and non- 
intersecting radiation. 

The present technique accommodates such intersection and non-intersection of 
radiation by generating a count map for each slice. In general, the count map is a plan 
25 view, on a pixel-by-pixel basis, wherein each pixel is assigned a value corresponding 

to the number of times it is traversed by radiation from the various projections (i.e. 
source locations). The count maps, then, represent the number of back-projection 
paths of the radiation through each pixel at each slice level. Fig. 4 represents an 
exemplary slice map for an elevated slice. The slice map 72 will generally include 
30 boundaries 74 defined by counting or computing the number of times each individual 

pixel is traversed by radiation in the various projections. The boundaries, then, define 

10 > 
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regions 76 in which pixels have the same count. Again, the values are assigned on a 
pixel-by-pixel basis for various pixels of the slice image as represented generally by 
reference numeral 78. 

5 Several aspects of the count maps are worth noting. First, it is generally 

anticipated that higher counts will be present in center regions of the map as compared 
to regions that are closer to the detector boundary, due to the fact that outlying source 
positions generally still traverse the central region of each slice. Similarly, it will be 
anticipated that lower slices would have larger regions of higher counts than upper 

10 slices due to the angles of the projections and the fact that the projections generally 

are directed toward the portion of the imaging plane in which the detector is located. 
As will be appreciated by those skilled in the art, this is true only for stationary 
detectors, to minimize the total "low count" area or volume, the detector could move 
such that the center slice has the biggest high-count region. Also, the particular 

15 contours and locations of boundaries 74 and of regions 76 may vary considerably 

depending upon the arrangement of the source positions giving rise to the various 
projections. 

In the present technique, the count map 72 for each slice is utilized to weight 
20 the contribution of pixels in each to the sensed X-ray attenuation. That is, regions of 

the count map having the highest value are considered to have fully contributed to the 
X-ray attenuation, where as counts that are lower are considered to contribute less, 
and are therefore weighted a lower amount in reconstruction. The weighting based 
upon the account map may be linear, that is, may vary proportionally with the count 
25 number, or it may also be non-linear. 

Further non-uniform weighting is provided by the present technique to account 
for a false assumption in conventional back-projection that results in contribution of 
regions entirely outside of the subject to sensed X-ray attenuation. Fig. 5 generally 
30 represents a slice weighting technique implemented by the present processing. As 

shown in Fig. 5, subject 18 is positioned immediately adjacent to detector 22, and 
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slices 66 are to be computed by back-projection. As will be appreciated by those 
skilled in the art, however, conventional techniques may result in ghost-like images 
being computed for regions far outside (e.g. to the left of) subject 18. The present 
technique employs slice weighting in which non-uniform weighting factors are 
5 computed, particularly outside the anticipated or actual bounds of the subject. In the 

illustration of Fig. 5, the height above the imaging plane adjacent to detector 22 is 
indicated along in axis 80, and a weighting factor value is indicated along a vertical 
axis 82. The weighting function is represented by a trace 84 which varies in value 
with the height above the detector. In the illustrated embodiment, an upper ramp 

10 weighting region 86 is defined which provides weighting values down to a null value 

at a desired distance from the upper limits of the subject. A similar lower ramp 
weighting region 88 is provided near the lower bound of the subject. A central 
weighting portion 90, which includes the slices of interest, has a uniform or generally 
uniform value of unity. It is anticipated that certain optimization of the non-uniform 

15 slice weighting values may be determined, and in the context of medical imaging, that 

the reduction from a value of unity to a null value may begin at, just below or just 
above the skin line of the subject. As discussed below, such non-uniform view 
weighting may be used with or separately from the count map weighting described 
above. 

20 

A further issue with a conventional tomosynthesis back-projection involves 
different system performance at different projections. Fig. 6 illustrates schematically 
this point. As shown in Fig. 6, for two different source locations, resulting in two 
different projections, an object 92 may appear substantially differently. In the 

25 example of Fig. 6, for example, a spherical object 92 may appear in a generally central 

projection 94 as generally circular, as indicated at reference numeral 96. As the 
source is moved to other positions to produce other projections, such as projection 98 
in Fig. 6, the same object may appear deformed or distorted, such as a spherical object 
projecting as a generally elliptical shape 100 in an outline projection as shown in Fig. 

30 6. In general, such distortions may be caused by a wide range of factors. These 

factors may be a function of the source position, the detector performance, particularly 
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near the edges of the detector, the geometry of the system, the distance traversed by X- 
rays in the various projections, beam scatter, beam hardening, and so forth. 
Conventional back-projection techniques, however, do not account for such system 
variation. All projections in conventional techniques are considered to be equally 
valid for all pixel locations when considered for computation of the appropriate pixel 
values in the various reconstructed slices. It should be noted that the projection 
weighting can also be used to manage artifacts. For example, an initial segmentation 
of the imaged object could be made (e.g., into chest wall and lung tissues), with the 
resulting regions receiving different weights (e.g., regions of interest, such as lung 
tissues receiving higher weights). This could reduce artifacts at the skin line in 
medical images. 

In the present technique, on the other hand, individual projection weighting 
maps are generated for each of the projections utilized to produce the reconstructed 
15 slice images. The projection maps may be generated as a calibration step, and stored 

for later use in reconstruction of the slice images. In general, it is contemplated that 
such projection maps will be generated for each projection utilized in the 
reconstruction technique, as indicated generally at reference numeral 102 in Fig. 6. 

The various non-uniform weighting approaches described herein may be used 
together or separately. In general, the resulting weighting will change the computed 
contribution of the various slice features, slices and projections to the computed pixel 
values for each re-projected slice image. The general relationship between the non- 
uniform weighting factors discussed above may be represented by the relationship 

I i = W s [I i W pi /W ci ]: 

where Ii is the value of each individual pixel in the reconstructed slice image, W s is 
the slice weighting factor for each individual slice, W P j is the projection map 
30 weighting factor for each individual pixel in each projection, and W C j is the count map 

weighting factor for each pixel in each individual slice. Note that (IiW P 0 has summed 
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up all contributions from different projections and is appropriately normalized. The 
use of the count map generally will cause a non-homogeneous look in the resulting 
images to be compensated for during back-projection. The use of slice weighting will 
cause out-of-plane artifacts due to limited projection angles to be suppressed for 
5 regions near or beyond the limits of the subjects of interest. Finally, the use of 

projection weighting will account for non-uniformities in the tomography system 
response. 

The concept of back-projection reconstruction with view weighting in 
accordance with the present techniques may be thought of as consisting of several 
steps. First, multiple digital tomosynthesis projections are acquired and stored, these 
projections are then accessed for processing. One, two and three-dimensional view 
weighting functions (masks or maps) for projection, slice and volumes of interest are 
then computed. The view weighting functions generally consist of spatially- varying 
values or weights. Examples provided in the present context include a count map 
where entries in the map represent the number of times that back-projection rays 
intersect pixels at various slice elevations. Functions, such as box functions, ramping 
functions, for smooth functions are applied for contributions of various slices to X-ray 
attenuation. Finally, overall system response is accommodated by production of 
projection maps used for weighting in the reconstruction process. 

While the invention may be susceptible to various modifications and 
alternative forms, specific embodiments have been shown by way of example in the 
drawings and have been described in detail herein. However, it should be understood 
25 that the invention is not intended to be limited to the particular forms disclosed. 

Rather, the invention is to cover all modifications, equivalents, and alternatives falling 
within the spirit and scope of the invention as defined by the following appended 
claims. 
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